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We determine the inclusive B — > X c lv branching fraction, B c e v , the CKM matrix element \V c b\, 
and other heavy-quark parameters from a simultaneous fit to moments of the hadronic-mass and 
lepton-energy distributions in semileptonic B-meson decays, measured as a function of the lower 
limit on the lepton energy, using data recorded with the BABAR detector. Using Heavy Quark 
Expansions (HQEs) to order 1/mf, we extract B cev = (10.61 ± Q.lQ exv ± 0.0Ghqe)% and |V C &| = 
(41.4 ± 0.4 exp ± OAhqe ± 0.6th) X 10~ 3 . The stated errors refer to the experimental, HQE, and 
additional theoretical uncertainties. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



The CKM matrix element V c b is one of the fundamen- 
tal parameters of the Standard Model and thus its pre- 
cise measurement with well understood uncertainties is 
important. At the parton level, the weak decay rate for 
b — > civ can be calculated accurately; it is proportional 
to \V c b\ 2 and depends on the quark masses, mi, and m c . 
To relate measurements of the semileptonic i?-meson de- 
cay rate to | V c b | , the parton- level calculations have to be 
corrected for effects of strong interactions. Heavy-Quark 
Expansions (HQEs) [l| have become a useful tool for cal- 
culating perturbative and non-perturbative QCD correc- 
tions [2| and for estimating their uncertainties. In the 
kinetic-mass scheme [3j for example, these expansions in 
l/rrib and a s (mb) (the strong coupling constant) to or- 
der 0(1/ ml) contain six parameters: the running kinetic 
masses of the b— and c— quarks, mt(/j) and m c (/x), and 
four non-perturbative parameters. We determine these 
parameters from a fit to the moments of the hadronic- 
mass and electron-energy distributions in semileptonic B 
decays to charm particles, B — > X c lv. This fit yields 
measurements of the inclusive branching fraction B c g v 
and of | V c b\, significantly improved compared to earlier 
BABAR measurements |4j. It also allows us to test the 
consistency of the data with the HQEs employed and to 
check for the possible impact of higher-order contribu- 
tions. Moment measurements and the extraction of |V c b| 
based on HQEs y| were first performed by the CLEO col- 
laboration [6j. More recently, global fits to a variety of 
moments were presented 00, 0, using HQEs in different 
mass schemes. 

This analysis makes use of moments measured by 
the BABAR collaboration [ToL Hl| . The moments are 
derived from the inclusive hadronic-mass (mx) and 
electron-energy (Eg) distributions in B — > X c lv de- 
cays, averaged over charged and neutral B mesons pro- 
duced at the T(4S) resonance. We have subtracted the 
charmless contributions based on the branching frac- 
tion S^=(0.22 ± 0.05)% [HI. All moments are mea- 
sured as functions of E cut , a lower limit on the lepton 
energy (for energy moments we only use electrons, for 
mass moments we also use muons). The moments are 
corrected for detector effects and QED radiation [l3| . 
The hadronic-mass distribution is measured in events 
tagged by the fully reconstructed hadronic decay of the 
second B meson. The hadronic-mass moments are de- 
fined as M*(E cut ) = {m x ) Ee >E cut with n =1,2,3,4. 



The electron-energy distribution is measured in events 
tagged by a high-momentum electron from the sec- 
ond B meson. We define the first energy moment as 
Mf(E cut ) ~ (Ee) Ei>E cut and the second and third mo- 
ments as Mi(E cut ) = {(Ei - Mi(E cut )) n ) Et>Ecut with 
n =2,3. In addition, we use the partial branching frac- 
tion M$(E cut ) = Jz2 a t *(dB ce „/dE e )dE e . All measured 
moments are shown in Fig. ^ 

In the kinetic-mass scheme the HQE to Q(l/mj) for 
the rate of B — > X c lv decays can be expressed as [l4| 



civ 



1927T 3 



V cb \ 2 (l + A ew )A pert (r,fi) x 



*o(r) 1 



,,2 .,2 , P'd+Pls 

2m 2 b 



-2(1 -r) 



.,2 j_ P'd+Pls 

4 Hq 



^+d(r)^ + 0(l/mi) 



(1) 



The leading non-perturbative effects arise at 0(1/ m 2 ) 
and are parameterized by /i^(/i) and ^(p), the expecta- 
tion values of the kinetic and chromomagnetic dimension- 
five operators. At C(l/m 3 ), two additional parameters 
enter, p^,(p,) and p\g(n), the expectation values of the 
Darwin (D) and spin-orbit (LS) dimension-six operators. 
These parameters depend on the scale fj, that separates 
short-distance from long-distance QCD effects; the cal- 
culations are performed for [i = 1 GeV [3j . Electroweak 
corrections are 1 + A ew = (\ + a/ir\xi(Mz /mb)) 2 = 1.014 
and perturbative QCD corrections are estimated to be 
A pert (r,n) = 0.91 ± 0.01 0. The ratio r = m 2 /m\ en- 
ters in the phase-space factor zo(r) = 1 — 8r + 8r 3 — r 4 — 
12r 2 bar and the function d(r) = 8 lnr + 34/3 - 32r/3 - 
8r 2 + 32r 3 /3 - 10r 4 /3. 

This analysis uses linearized expressions for the 
HQEs 15]. Specifically, the dependence of \V c b\ on the 
true values of heavy-quark parameters, expanded around 
a priori estimates of these parameters, is 



B 



civ 



1.55 ps 



0.0417 y (0.1032) t b 

x[l + 0.30(a s (m fc ) -0.22)] 

x [1 - 0.66(m b - 4.60) + 0.39(m c - 1.15) 

+0.013(/4 - 0.40) + 0.090!, - 0.20) 

+0.05(/4 - 0.35) - 0.01(pi s + 0.15)]. (2) 
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FIG. 1: The measured hadronic-mass (a-d) and electron-energy (e-h) moments as a function of the cut-off energy, E cut , 
compared with the result of the simultaneous fit (line), with the theoretical uncertainties |15f indicated as shaded bands. The 
solid data points mark the measurements included in the fit. The vertical bars indicate the experimental errors; in some cases 
they are comparable in size to the data points. Moment measurements for different E cu t are highly correlated. 



Here m^, and all other parameters of the expansion are 
in GeV^ ; tb refers to the average lifetime of B mesons 
produced at the T(AS). We use tb — fo T o + (1 — fo) T ± = 
(1.608 ± 0.012) ps, taking into account the lifetimes [lfj 
of neutral and charged B mesons, tq and t±, and their 
relative production rates (defined in terms of /o = 0.488± 
0.013 H3, the fraction of B°B° pairs). 

HQEs in terms of the same heavy-quark parameters 
are available for the hadronic-mass and electron-energy 
moments. The dependence on the heavy-quark param- 
eters has again been linearized u sing the same a priori 
estimates of the parameters [14], [l5|. We have verified 
that the differences between the linearized expressions 
and the full theoretical calculation are very small in all 
cases. We use these linear equations to determine the 
unknown heavy-quark parameters, the total branching 
fraction B c i Vl and \V c b\ from a simultaneous % 2 fit to the 
measured moments and the partial branching fraction, 
all as a function of the cut-off lepton energy, E ' cut . 

In total, we have measured four hadronic-mass mo- 
ments for each of seven different values of E cut , ranging 
from 0.9 to 1.5 GeV, and three electron-energy moments 
plus the partial branching fraction at five values of E cu t, 
ranging from 0.6 to 1.5 GeV Since many of these 

individual moments are highly correlated we select for 
the fitting procedure a set of moments for which the cor- 
relation coefficients do not exceed 95%. Thus we only 
use half of the 28 mass moments, and retain 13 of the 20 
energy moments. 

The global fit takes into account the statistical and 



systematic errors and correlations of the individual mea- 
surements, as well as the uncertainties of the expres- 
sions for the individual moments. We assess the un- 
certainty of the calculated moments by varying, as sug- 
gested in [l^, in the linearized expressions (given for 
| Vet | in Eq. EJ) the a priori estimates for /z£ and p? G 
by ±20% and for p 3 LS and p 3 D by ±30%. We assume 
that for a given moment these changes are fully corre- 
lated for all values of E cut , but uncorrelated for different 
moments. The resulting fit, shown in Fig. ^ describes 
the data well with \ 2 = 15-0 f° r 20 degrees of free- 
dom. Tables [I] and [H] list the fitted parameters, their 
errors and correlations. As expected, and m c are 
highly correlated. For the mass difference, we obtain 
m b ~m c = (3.436 ±0.025 exp ±0.018 i? Q B ±0.010 a ,) GeV. 

Beyond the uncertainties that are included in the 
fit, the limited knowledge of the expression for the de- 
cay rate, including various perturbative corrections and 
higher-order non-perturbative corrections, introduces an 
error in \V c h\, assessed to be 1.5% ^bi|. By compari- 
son, the impact of the uncertainty in a s is estimated to 
be relatively small. For M^(E cut ) moments, perturba- 
tive corrections of order a s are included with a s {m,b) = 
0.22±0.04, whereas for M^{E cut ) moments, they are cal- 
culated only to 0[a 8 ) with a s (rrib) = 0.3 ±0.1. We esti- 
mate the error on the perturbative corrections by varying 
a s within the stated errors. The choice of the scale \x is 
estimated to have a very small impact on \V c b\ and the 
branching fraction [l4T |. 

A series of tests has been performed to verify that the 
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TABLE I: Fit results and error contributions from the moment measurements, approximations to the HQEs, and additional 
theoretical uncertainties from a s terms and other perturbative and non-perturbative terms contributing to T^. 





\v cb \(w~ 3 ) 


m b (GeV) 


m c (GeV) /4 


(GeV*) p A D (GeV a ) 


P-g (GeV*) pi s (GeV a ) 


Bah, (%) 


Results 


41.390 


4.611 


1.175 


0.447 


0.195 


0.267 


-0.085 


10.611 




0.437 


0.052 


0.072 


0.035 


0.023 


0.055 


0.038 


0.163 


&HQE 


0.398 


0.041 


0.056 


0.038 


0.018 


0.033 


0.072 


0.063 




0.150 


0.015 


0.015 


0.010 


0.004 


0.018 


0.010 


0.000 


Sr 


0.620 
















Slot 


0.870 


0.068 


0.092 


0.053 


0.029 


0.067 


0.082 


0.175 



TABLE II: Correlation coefficients for the fit parameters. 
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in i, m c fj%~ pi, Pg Pls B c t v 



Pa 



Pls 



\Vcb\ 

m b 
m c 

A 

I'G 

Pls 

Bclv 



1.00 



-0.49 
1.00 



-0.36 
0.97 
1.00 



0.56 
-0.40 
-0.38 

1.00 



0.35 
-0.13 
-0.13 
0.82 
1.00 



-0.37 
0.16 

-0.04 
0.08 
0.08 
1.00 



0.64 
-0.63 
-0.50 

0.46 

0.23 
-0.43 

1.00 



0.61 
0.23 
0.29 
0.16 
0.12 
-0.04 
0.09 
1.00 



fit results are unbiased. Specifically, we enlarged and re- 
duced the estimated theoretical uncertainties by a factor 
of two and verified that the changes in the fitted param- 
eters were small compared to the errors of the standard 
fit. We have also checked that the choice of the set of 
moments that are used in the fit does not significantly af- 
fect the result. In particular, an energy cut-off above 1.2 
GeV might have introduced larger theoretical uncertain- 
ties and a potential bias, and moments for lower values 
of E cut might have been affected by higher backgrounds. 
We found no evidence for any such effects. 

The fit results are fully compatible with independent 
estimates |l5( of /ig=(0.35 ± 0.07) GeV 2 , based on the 
B* - B mass splitting, and of p\ s =(-Q.lb ± 0.10) GeV 3 , 
from heavy-quark sum rules |18j . 

Figure shows the A\ 2 = 1 ellipses for \V c b\ ver- 
sus mj, and /i 2 , for a fit to all moments and separate 
fits to the electron-energy moments and the hadronic- 
mass moments, but including the partial branching frac- 
tions in both. The lepton-energy and hadronic-mass mo- 
ments have slightly different sensitivity to the fit pa- 
rameters, but the results for the separate fits, \V c b\ = 
(41.4 ±0.7) x lO- 3 and \V cb \ = (41.6 ± 0.8) x 10" 3 are 
fully compatible with each other and with the global fit 
to all moments. Changes in the other fit parameters are 
also consistent within the stated errors. Since the expan- 
sions for the two sets of moments are sensitive to different 
theoretical uncertainties and assumptions, in particular 
the differences in the treatment of the perturbative cor- 
rections, the observed consistency of the separate fits in- 
dicates that such differences are small compared with the 
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FIG. 2: Fit results (crosses) with contours corresponding to 
A% 2 = 1 for two pairs of the eight free parameters a) nib and 
b) p% versus \ V c b\, separately for fits using the hadronic-mass, 
the electron-energy, and all moments. 



experimental and assumed theoretical uncertainties. 

In conclusion, we have extracted \ V c b\, the semileptonic 
branching fraction, and the heavy-quark masses, 



(41.4 ±0.4 exp ± QA HQE ±0.6*) 

To, 



x 10" 



\V cb \ = 

B ceu = (10.61 ±0.16 e:cp ±0.06 ff Q £ ) 

m fc (lGeV) = (4.61 ± 0.05 exp ± 0M H qe ± 0.02 th ) GeV, 

m c (lGeV) = (1.18 ± 0.07 exp ± 0.06hqe ± 0.02 t/l ) GeV, 

as well as the non-perturbative parameters in the kinetic- 
mass scheme up to order (l/m\) (see Table QJ. The total 
semileptonic branching fraction is B cev + B ueL , = (10.83 ± 
0.16 exp ± 0.06hqe)%- The errors refer to contributions 
from the experimental errors on the moment measure- 
ments, the HQE uncertainties included in the fit, and 

additional theoretical uncertainties, 8 t h 

derived from Refs. 0G3. 

Based on a large set of hadronic-mass and electron- 
energy moments and a consistent set of HQE calcula- 
tions, we have also been able to assess the uncertainties 
in the 0(1/ m 3 ) terms from the data without constraints 
to any a priori values. The fitted values of the parame- 
ters are consistent with theoretical estimates HQ. The 
uncertainties on the quark masses are much smaller than 
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those of previous measurements |16]. It would be inter- 
esting to compare the results of this analysis with fits 
based on recent calculations performed in the IS mass 
scheme [il^. 

The result on \V c b\ is in agreement with previous mea- 
surements using HQEs, either for a different mass scheme 
and with fixed terms of C(l/m^) ||, or for the kinetic- 
mass scheme, but with external constraints on almost all 
HQE parameters @ , as well as with an analysis combin- 
ing both of these measurements Q • 
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